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Abstract 



i-C ' An analysis including most recent Belle data on A (3872) is per- 

formed, using coupled channel Flatte formula. A third sheet pole close 
to but below D^D*^ threshold is found, besides the bound state/virtual 
(— I I state pole discussed in previous literature. The co-existence of two poles 

near the D^D*^ threshold indicates that the A (3872) may be of ordinary 
^ , cc 2^Pi state origin, distorted by strong coupled channel effects. The 

^> ' latter manifests itself as a molecular bound state (or a virtual state). 
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§ ; 1 Introduction 

In year 2003 the Belle collaboration found a very narrow [Fx < 2.3 MeV) 
^ I resonance structure named X(3872) in the J/^tttt invariant mass spectrum, 

I in the J/\E' tt+tt^ process [1]. The branching ratio Br(i?+ 

K+X) Br(A n+TT^J/^) is updated to be = (7-10) x lO'^ both by BaBar [2] 
and by Belle [3]. Moreover, Belle also observed A(3872) in the decay and 
found that the rate is comparable with that of the charged channel [3]. Most 
recently, the CDF Collaboration reported a new measurement on the mass 
parameter in the J/'^ir'^ir" channel [1], 

Mx = 3871.61 ±0.16 ±0.19 MeV . (1) 

Replacing the old CDF measurement by the new one results in a world average 
of Mx = 3871.51 ± 0.22MeV, which is very close to the D°D*° threshold 
M^o^.o = 3871.81 ±0.36 [0]. 
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The other decay modes of X(3872) include J/^'7r+7r^7r° [li], [7] and 

\l/'7 [7] with relative rates 

Br{X vr+vr- J/*) ' ^ ^ 



Br(X ^ 7J/^) 
Br(X ^ n+TT-J/^) 

Br(X ^ 7$') 
Br(X ^ TT+vr- J/^) 



0.33 ±0.12, (3) 
1.1 ±0.4. (4) 



The dipion mass spectrum in the J/\E'7r"'"7r~ mode shows that they come from 
the p resonance [1] and the Stt in the J/'^ir'^n'ir^ mode come from the u 
resonance [6]. Thus, the rario i? ^ 1 in (2) indicates that there should be large 
isospin violation in the decays of X(3872). 

In year 2006 the Belle collaboration studied the D^D^t^^K^ decay 

process and found the enhancement of the D^D^n^ signal just above the D^D*^ 
threshold [8], the resonance is peaked at 

Mx = 3875.2 ± 0.7^?:^ ± 0.8 MeV, (5) 

roughly 3.6 MeV higher than the value in (1). The corresponding branching 
ratio at the D^D^-k^ peak is [8], 

Br(5+ ^ K+D''D%'') = (1.02 ± Mtoll) x 10"' . (6) 

The different peak locations of X(3872) in the D^D^tx^ (D^D*^) and J/^Tr+Tr" 
channels are reconfirmed by latter BaBar experiments [9]. In 2008, a new 
analysis to the Belle data in the {D*° D^tt^ and D*^ D^-f) 

channel is given [10], and the new determination of the peak is 2.6MeV lower 
than the previously reported by Belle [10], 

Mx = 3872.6+°:^ ± 0.4MeV . (7) 

The difference comes from the inclusion of new data {D* D'-f), more so- 
phisticated fit (unbinned fit with mass dependent resolution), and improved 
Breit-Wigner formula (the Flatte formula). The central value as given by 
Eq. (7) is, however, still about IMeV above than the value measured by CDF 
Collaboration []]. Meanwhile in Ref. [10] a renewed determination of the fol- 
lowing branching ratio is given, 

Br(5+ ^ K+X{D*^D^)) = (0.73 ± 0.17 ± 0.13) x 10"^ . (8) 

The X(3872) is naturally interpreted as a C = ± molecule of D^D*^ in 
s-wave [11, 12] since its mass is very close to the D^D*^ threshold and the 
quantum number J^'-' = 1+"*" is favored by the experimental analysis [13]. 
It also predicted the J/'^uj mode with similar rate as J/'^p [12]. However, 
the large production rates of X(3872) in B-factories and at Tevatron favor a 
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conventional charmonium assignment [14, 15] (say, Xci) rather than a loosely 
bound state of D'^D*^. Furthermore, the large decay rate of X — >■ ^^'7 in (4) 
also strongly disfavors the molecular assignment since it is very difficult for the 
transition of a molecular to \E^' through the quark annihilation mechanism [12]. 
The large isospin violation indicated by (2) can also be explained quite well in 
the charmonium model [16]. Hence it seems that we are facing a dilemma in 
recognizing X(3872). 

To further clarify the identity of X(3872), one needs to look deeper into 
the pole structures of the scattering amplitude involving X(3872). For a dy- 
namical molecule of D^D*^, there is only one pole near the threshold, and the 
requirement of two nearby poles to describe the X(3872) will generally imply 
that it is a cc state near the threshold [17]. The line shapes of — >■ XK^ in 
the J/\Ef-7r+7r~ and D^D^tt^ / D^D*^ modes and the corresponding pole struc- 
tures have been studied by two groups [18,19] independently. Both fits give an 
one-pole structure, although one fit [IS] favors a virtual state and the other [19] 
favors the loosely bound state. Since more data are available after these two 
fits, it deserves a careful reanalysis to the data of X(3872). In this paper we 
devote to the study of this problem. In Sec. 2, we firstly describe the method 
we use for the analysis, we also describe how we make the fit from various 
experimental data. Sec. 3 is for the discussions and conclusions. The final 
result of this analysis presents a unified picture in understanding the dual faces 
of X(3872): A third sheet pole close to but below D^D*^ threshold is found, 
besides the bound state/virtual state pole discussed in previous literature. The 
co-existence of two poles near the D'^D*^ threshold indicates that the X(3872) 
may be of ordinary cc 2^Pi state origin, distorted by strong coupled channel 
effects. The latter manifests itself as a molecular bound state (or a virtual 
state). 

2 Coupled channel description of the X(3872) 
resonance 

Notice that X(3872) associates with nearby different cuts, hence a coupled 
channel analysis is needed in order to take care of the complicated singularity 
structure. This has already been emphasized in Refs. [18,19,20,21]. Hanhart 
et al. gave a very interesting explanation to the X(3872) peak as a virtual 
state [18]. Their conclusion relied on of course the experimental data available, 
and especially on the two peak structure in different channels. The latter plays 
a crucial role in getting such a conclusion. In the analysis of Hanhart et al., the 
effect of energy resolution is not considered. Since the two peaks are not too far 
from each other and the difference between them is comparable in magnitude 
to the energy resolution parameter, one worries about that the negligence of 
energy resolution effect may distort their conclusion. For reasons mentioned 
previously a new analysis on this subject is necessary. We proceed with data 
prsently available [2,3,9,10] to reanalyze the X(3872) state, with the energy 



3 



resolution effect talcen into account. On the theory side the method we use is 
essentially the same as that of Ref. [18]. 

For describing the chain decays with X(3872) involved as intermediate 
state, we parameterize the inverse of the propagator of X(3872) as 

D{E) = E-Ef+ '-{g,h + 92k2 + T{E) + , (9) 



where Ej = Mx - Mno - M^.o; h = yfZjME, k2 = y2/i2(^ - 5) and 6 = 
Md+ + Mo*- — Mjjo — M^.o, /ii and fi2 are the reduced masses of D^D*^ 
and D~^D*~, respectively. Isospin symmetry requires gi ~ g2. ^{E) includes 
channels J/^tt+tt^ (through J/'^ p), J/'^tv^ti^ti^ (through J/'^uj) channels. 
Different from Ref. [18] here we add a constant width Tc to simulate every other 
channels, including radiative decays and light hadron decays. From Eq. (4) we 
know that this term is certainly non- negligible as comparing with J/'^'k^t[~ 
decay, not to mention the to be observed light hadronic decays. 

For simplicity, we describe the p and uo resonances in the final states by 
their Breit-Wigner distribution functions, then one has 

.Mx-rnjf^ (im k{m)Tp 



^w+w-J/'fiE) = fp 

J2 



2m,r 27r (m — mp)2 + r2/4 



fMx-mj/^ dm k{m)T^ 
W.o,/^(E) = UJ^^^ 27(m-mJ2 + r2/4' ^^^^ 

where fp and are the X couplings to J/'^p and J /'^u respectively, Mx = 
E + M£,o + Mj^*o is the (off-shell) center of mass energy of the X particle and 



k{m) 



{Ml - (m + mj,^Y){Ml -{m-mj,^y) 



\ AM'- 



Let B = Br{B XK), recalling that Br{D*^ D°7r°) = 61.9±2.9%, [22] 
repeatedly using the chain decay formulae leads to, 

dE = °-^'^2^ \D{EW ' 

dBT[B Kn+Ti-J/^ _ ^1 T^+^-j/^{E) 



dE 271 \D{E)\ 
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dBilB ^ JiTT+TT-Tr" J/^] ^ 1 T^+^-^oj/^jE) 

dE 271 \D{E)\^ ■ ^ ' 

In the fit to X —s> D*^D^ data [9], since all decay modes of D*^ are considered 
there, we drop the factor 0.62 in the first formula of the above equation. 

One also has to consider the background contributions. In all the fit to the 
data, we assume there is no interference between data and background. This 
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is in coincidence with experimental analyses. In DDir channel we assume the 
background contribution is proportional to EoDn, hence 

For the D*^D^ final state we assume background contribution is proportional to 
the phase space of D^D*^, ki. In the J/\E'7r+7r~ case, we assume the background 
is a constant. Herewith we often use overlined branching ratios to represent 
the signal plus background contributions: 

+ (14) 

dE dE ' ^ ^ 

The ratio R defined in Eq. (2) has to be put into the fitting program as a 
constraint. Throughout this paper the ratio R = ^- "^^^ 

formula used to estimate R is the same as that adopted by Hanhart et al. [18], 
but in here we constrain the value R by using the penalty function method, 
which is simple and effective. That is we effectively add a term to the total 
X^'- X/? = 10 X |i? — lp/0.4^. The factor 10 is an arbitrarily chosen penalty 
factor which is enough to make the ratio R being almost exactly unity. We 
notice that the ratio R measured by experiments contains a large error bar as 
shown in Eq. (2). We will therefore also pay some attention in the numerical 
fit to different value of R, in next section. 

3 The data fitting program and the fit results 

3.1 Data samples and the energy resolution parameters 

As stated earlier we use 4 sets of data: 

1: The X D*^D^ mode by BaBar [9], where D*^ is reconstructed both 
from D^TT^ and D^'j mode. There are 12 data points in the fit region 
from D*^D^ threshold up to 3.895GeV. The background contribution 
starts from D*^D^ threshold, the same as that adopted in Ref. [9]. The 
corresponding number of events distribution is, 

N^"^*' - 2fMeVl ^^-^ rfB?[i? KD^D*^] 

2: The X D^D^ir^ data from Belle Collaboration [S] is replaced by the 
upgraded one from — > XK^ [10]. We fit the data in the energy re- 
gion from D^D^Ti^ threshold to 3.91257GeV, there are totally 119 events 
collected from decays. 

3: Data of J/^tt^tt" from BaBar [2]. We only use the charge mode {B^ 

X(3872)i^:+) data, since the error bar of the neutral mode (5° ^ X(3872)i^°) 
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data are much larger. There are 11 data points in the fit region 3.84 < 
Mx < 3.89GeV and 



NsaBar = 5[MeV] ^ . (16) 

4: Data of J/\l/7r+7r" from most recent Belle experiments [3]. We fit the 
data sample in the energy region from 3.84135GeV to 3.90173GeV, with 
totally 398 events. 

With the unbinned data sets from Belle Collaboration on both X — > 
D^D^n^ from XK^ and X J/'^t\:^'k~ decay, we make a combined 

fit of likelihood method and method in the following way: 

Xlff = -2 E log + E X? + , (17) 

where i = 2, 4; j = 1, 3. The background contributions to the two data samples 
2 and 4 are treated similarly as those discussed previously. The PDF used in 
the likelihood fit is written as, 



+ h.g. 
JdE \'^ + b.g. 



dE 



Because the peaks in different channels are rather close to each other, one 
needs to take energy resolution effect into account. 



Br(E)=jiE^Br(E,)^=— (19) 



In general, the energy resolution parameter o" is a function of Er^, the original 
energy of incoming particles. For J/\E^7r+7r" channel at Belle: (t{Ex) = 3MeV . 
For D*°D^ at Belle: a{E^) ~ 0.176^^^ - Md*odo ■ [10] We assume that the 
BaBar detector maintains the same energy resolution parameters. 



3.2 Pole locations determined from combined data fit 

Experiments indicate B to be about a few times 10~^. The value of B is about 
2-4 X 10~^ in the charmonium model [14], while in the molecular model, it is in 
general not larger than 1 x 10~^ [19,23]. Therefore in the following analyses, 
we often fix i3 at a few times 10~^, though it is noticed that the fit program 
prefers a larger value of i3 ~ 2 x 10^^ with large error bars. 

We have stressed in section 2 that we add a constant width term in the 
Flatte propagator, which corresponds to modes rather than the near threshold 
ones (J/\Efp, J/'^u, DDn). These modes include both the observed ones, 
such as '^^'^'y [7], and the hidden ones. In the charmonium model, the most 
important hidden decay mode of X(3872) as Xci(2P) is the inclusive light 
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hadronic decay, and the partial width is of 0{1) MeV [16]. However, for the 
pure D^D*^ molecule, it is difficult to annihilate the charm quark pair into 
light hadrons. Therefore, the most important hidden modes of X(3872) in the 
molecular model may be the hadronic transitions to Xcj(l-P); such as Xco'^^ 
and XciTTTT, while the widths of them are expected to be smaller than that of 
J/\i/7r+7r~ [24]. Thus, the term Fc can provide important information on the 
X(3872). 

Poles on different sheets are searched for using results of fit parameters. 
The naming scheme of Riemann sheets is given in table 1. In table 2 - 5 we 





II 


III 


IV 








+ 




+ 







Table 1: Naming scheme of Riemann sheets. 



list several fit results with different choices of i3 ~ a few x 10~^. The error 
of Fc is given while others are not listed. For comparison we also list the fit 
results by setting Fc = 0. Notice that in tables 2-5 the parameter gx relates 
to parameter gi in Eq. (9) as g, = ^^^^jl^^^^^2 . 



B = 2x 10-^ 


gxiGeV) 


Ef{MeV) 


fp X 103 


X 102 


Fe(MeV) 


Xeff = 4092 


4.16 


-6.79 


2.10 


1.45 


1.78 ± 1.66 


X^ff = 4093 


4.40 


-6.40 


0.44 


0.32 





Table 2: Pole locations: Ej/^ = M - zF/2 = -4.72 - 1.51zMeV, Ej/ = 
M-zF/2 = -0.20-0.38zMeV (with FJ; Ej/^ = M-iV/2 = -3.72-.08iMeV, 
Ej^^ = M - iV/2 = -0.02 - O.OlzMeV (w/o FJ. 



i3 = 3 X 10-^ 


gxiGeV) 


Ef{MeV) 


fp X 103 


U X 102 


rc(MeV) 


= 4090 


4.20 


-6.89 


1.46 


1.01 


2.02 ± 1.61 


= 4092 


5.57 


-10.3 


0.74 


0.53 





Table 3: Pole positions: Ej/^ = M - iT/2 = -4.82 - 1.582MeV, Ej/ = 
M-iF/2 = -0.20-0.40iMeV (with F,); Ej^^^ = M-iV/2 = -7.66-0. 122MeV, 
Ejl = M - iT/2 = -0.02 - O.OliMeV (w/o FJ 



By examining the numerical results as given in tables 2-5 we have the 
following observations: 

1. A third sheet pole is always found. When B gets large (~ 1 x 10"^)^ 
the pole locates far away below the D^D*^ threshold. In such a case it 
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B = 5x 10-4 


gx{GeV) 


EfiMeV) 


fp X 103 


f. X 102 


re(MeV) 


= 4088 


5.41 


-11.1 


2.05 


1.39 


3.23 ±2.54 


= 4091 


7.45 


-18.3 


1.67 


1.18 





Table 4: Pole positions: Ej/ = -0.13 - 0.39iMeV, Ej/^ = -9.20 - 2.54iMeV 
(with rj; Ej/ = +0.04 - 0.08iMeV, E^^^^ = -18-9 - 2.82iMeV (w/o TJ. 



B=lx 10-3 


gx{GeV) 


EfiMeV) 


fp X 103 


X 102 




= 4086 


6.28 


-15.4 


2.06 


1.38 


5.67 ± 1.04 


= 4090 


10.4 


-36.3 


3.49 


2.44 





Table 5: E^^ = M - iV/2 = -0.16 - 0.58iMeV, = M-iT/2 = -14.70 - 
4.30iMeV (with TJ; Ej/ = M-ir/2 = -0.44-0. 12iMeV, Ej^^^ = M-iV/2 = 
-35.5 - 1.40iMeV (w/o 



is understood that the whole data may well be fitted by a parametriza- 
tion with a single pole. In this sense, the X(3872) may be regarded as 
'dynamically generated'. However, for more reasonable choices of (i.e., 
smaller) B, the third sheet pole is rather close to the threshold and is 
certainly physically relevant. 

2. When B is small, the fit predicts a value of compatible with the 
quark model prediction on 2^ Pi state light hadronic decay width, i.e., 
~ IMeV [16]. 

3. The location of the pole near D*^D^ threshold is not stable in the sense 
that it may either locate on sheet II or sheet IV. The former corresponds 
to a D*^D^ molecule, the latter corresponds to a virtual state. The 
current analysis is not able to make a definite conculsion on the two 
scenario, though the former is more preferable. 

Besides above observations, in the fit when Tc is set to zero, we also confirm 
the approximate scaling law among different parameters [18]. However, when 
the constant width is added, the approximate scaling law no longer exists. 
Because of the approximate scaling law, the authors of Ref. [18] fix one of the 
parameters [g). Their choice is similar to the situation of table 5, corresponding 
to a large B. Hence it explains why in the analysis of Ref. [18] the discussions 
on the third sheet pole is missed, since the latter is quite distant away from 
the physical region under concern. However, a choice of i3 ~ 1 x 10-3 seems to 
be too large to be realistic. A third sheet pole is of typical resonance behavior 
and can be identified as the missing 2^Pi cc state. The puzzle remained here 
is why the third sheet pole locates below the D^D*^ threshold. A pole with 
such a behavior is sometimes called a "crazy resonance" . [25] 
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From tables 2-5 we notice that the location of the nearby pole is not very 
stable numerically, though it seems to prefer to locate on the second sheet 
(hence a molecule). The second sheet pole may however shift above the D^D*^ 
threshold, or even switch to sheet IV. Therefore we hesitate to make any 
definite conclusion on the location of this pole. The only solid statement 
that can be drawn from above numerical analysis is that the third sheet pole 
moves towards the D^D*^ threshold and hence becomes non-negligible when 
parameter B is within a few times 10~^. 





Figure 1: Fit results with B = 3 x 10~^, Fc free. In order to compare with 
data samples II and IV, we set both bin size to be 3.5MeV, and for the former 
we give the histogram. 

In above an effective minimization procedure with mixed function and 
likelyhood function is being used. To check weather the qualitative picture 
revealed depends on the particular choice of Eq. (17) or not, we also tested 
the binned data fit by taking 1 bin=3.5MeV and fit to the same energy region. 
It is found that the major conclusion of our qualitative result is unchanged - 
that a twin-pole structure is needed when B is small. Taking B = 3 x 10~^ 
for example, the pole locations are found to be: = M — iV /2 = —3.84 — 

l.THMeV, = M - iT/2 = -0.10 - 0.43iMeV (with Fe), to be compared 
with the results of table 3. 
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The influence of different choices of the value of R defined by Eq. (2) is 
also tested. Setting R = 0.5 and 1.5 for example, it is found that, in all two 
cases, a third sheet pole a few MeV below DD* threshold is found, except the 
second sheet pole very close to the threshold. Therefore the variation of ratio 
R does not distort the qualitative picture obtained in our numerical analysis. 

4 Discussions and conclusions 

It is actually not surprising that our analysis finds two poles - the occurrence 
of two poles is an intrinsic character of the coupled channel Flatte propagator. 
The importance of the current analysis is that, as it points out, for reason- 
ably chosen value of B, the third sheet pole locates quite close to the D^D*^ 
threshold and hence be physically relevant, except for the sheet II (or sheet IV) 
pole emphasized in previous literature. This picture is found to be unaltered 
when varying the fitting method and the value of R. The conclusion certainly 
depends on the simultaneous fit to experimental data in two channels. The 
statistics of data set I and II are not as good as the J/\l/7r"'"7r~ data, hence fu- 
ture improvement on experimental data in D^D^* and D^D^tt^ channels would 
certainly be helpful in clarifying the issue further. 

The two pole structure of the X(3872) state as revealed in this study is 
important, as we believe, in understanding correctly the nature of the X(3872) 
resonance. In this aspect, it can be helpful to learn some lessons from previous 
studies on the /o(980) resonance. Generally one pole structure was considered 
as crucial evidence in supporting the molecule identification of the /o(980) 
state in the literature. On the other side, the existence of two poles close 
to the threshold was often interpreted as an evidence against the molecular 
state origin of the /o(980) resonance [17]. Early studies of the /o(980) tend 
to identify it as having only one pole near the KK threshold, and hence a 
molecular state. [26] It was found later that the tttt, KK scattering data are 
much better described by allowing two poles near the KK threshold. [27, 28] 
In this picture, the third sheet pole may contain a large qq component, that its 
position close to the KK threshold is due to the attractive interaction in the 
KK channel. The sheet II pole is mainly of KK molecule nature. The X(3872) 
situation should be rather similar to the /o(980) case, except that in here the 
pole locations are distorted more severely by coupled channel effects. What 
we would like to stress here is that the two pole structure of the state X(3872) 
may reveal its dual faces: it is of cc origin due to the existence of the sheet 
III pole, but the coupled channel effect also manifests itself by presenting an 
additional pole, close to the D^D*^ threshold. The latter can also be explained 
as molecular bound state/virtual state. 

It should be stressed that, a pure molecular assignment of X(3872) en- 
counter a difficulty: the favored value of B lead to the width Fc to be roughly 
of 0{l)MeV. The pure molecular assignment of X(3872), however, would pre- 
dict a much smaller value of Fc as mentioned earlier. Thus, our analysis 
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supports that X(3872) is a mixing state of x'd and D^D*^ components [1 1,15]. 
A nearby x'ci below D*^D^ threshold is actually reported by quenched lattice 
QCD calculation. [29] The gap between the mass of x'ci ^be quark model [30] 
and the experimental one in Eq. (1) can be reduced when coupled channel ef- 
fect is taken into account. [31] Here it is worth emphasizing that the shift in 
the mass of a 'pure' x'ci ^be attraction of the D^D*^ threshold, not 

because of its mixing with other cc state. 

To conclude, the analysis given in this paper suggests the following scenario 
for X(3872): Firstly, there exists a sheet 11 (or sheet IV) pole very close to the 
D*ojjo threshold, this confirms previous results in the literature. Secondly, a 
fit to the data with a reasonable choice of B parameter requires the existence 
of a third sheet pole, but below D*^D^ threshold - this observation is new. 
With the uncovering of the existence of two poles a clear understanding on the 
ambiversion of X(3872) emerges - that it can be identified as a 2^Pi cc state 
strongly distorted by coupled channel effects. 
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